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Abstract. Quantum annealing is a generic solver of classical optimization problems that 
makes full use of quantum fluctuations. We consider work statistics given by a repetition of 
quantum annealing processes by employing the Jarzynski equality proposed in nonequilibrium 
statistical physics. In particular, we analyze a nonequilibrium average of the exponentiated work 
performed by a transverse field. A special symmetry, gauge symmetry, leads to a non-trivial 
relationship between quantum annealing toward different targets in the theory of spin glasses. 
We believe that our results will be a step toward an alternative realization of efScient quantum 
^ , computation as well as our better understanding of nonequilibrium behavior of systems under 

^ ' quantum control. 
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in 

fSj ■ 1. Introduction 

, We should often solve difficult problems in a reasonable time, which is closely related to the 

efficiency in cost and time in industry and distribution systems. A well-known instance is 



(N . to find the minimum path for a person to visit all the cities, called the traveling salesman 

' problem, which is a typical optimization problem [HE]- The goal of quantum computation is 

to solve problems hard to solve quickly on classical computer in a moderate time by use of 
superpositions, tunneling effects, and entanglement of quantum nature. Quantum annealing is 
- j , a method to realize quantum computation O |3l El [6]. Conventional quantum computation 

' uses successive, discrete unitary operations to enhance the probability of the desired state 

, representing the answer of a hard optimization problem. In contrast, quantum annealing 

uses continuous dynamics represented by a time-dependent Schrodinger equation. Dynamical 
behavior of quantum systems during such processes is still poorly understood, and an adiabatic 
control of quantum dynamics is often employed, called the quantum adiabatic computation [7J. 
Such a quantum adiabatic computation, however, has a bottleneck for a certain type of the 
optimization problems [8', '9] . 

A recent study of quantum annealing indicates a possible deviation from adiabatic controls 
|10j . In the present paper, we proceed along this line and analyze the work performed by 
quantum fluctuations in a non-adiabatic control of system parameter. The analysis makes 
use of several remarkable properties found in statistical physics. One is an exact relationship 
between nonequilibrium processes and equilibrium states at the beginning and end, the Jarzynski 
equality |1H [T2] . Another is a special symmetry, called the gauge symmetry, found in spin glasses 



\13\ I14j. As detailed below, a combination of these theories enables us to obtain exact relations 
between equilibrium and nonequilibrium quantities and to evaluate several quantities related to 
dynamical processes in quantum systems [15]. 

After an introduction of the Jarzynski equality, we develop our analyses by use of gauge 
symmetry to show the derivation of an exact identity related with the performed work during 
quantum annealing. 

2. Quantum annealing and spin glasses 

We consider the following time-dependent Hamiltonian in quantum annealing 

H{t) = -f{t)jY^njata^ - r (1 - fit)) Y: erf, (1) 

where f{t) = t/T, and t runs continuously from to T. The initial Hamiltonian H{0) is 
composed only of the transverse-field term —TJ2i '^f ■ The final Hamiltonian H(T) is chosen to 
describe the optimization problem to be solved. We take a simple spin glass Hamiltonian as 
H(T), the lb J Ising model, on an arbitrary lattice throughout this paper for convenience. The 
sign of the interaction Tjj follows the distribution function 

P(r,) = pHI - r,) + (1 - mi + r,) = (2) 

where p is the concentration of the ferromagnetic interaction J > 0, and exp(— 2/?pJ) = {l—p)/p. 
A sufficiently slow decrease of the strength of the transverse field changes the trivial initial state 
to the nontrivial ground state of the target Hamiltonian H{T). This is the idea of quantum 
adiabatic computation realized in T ^ oo. Instead of the adiabatic control, in the present study, 
let us consider the repetition of fast quantum annealing (T small) starting from an equilibrium 
ensemble. We may not be able to always reach the desired state since the system does not 
trace the instantaneous ground state when the adiabatic condition is not satisfied. Therefore we 
need to repeat the process to hit the correct ground state in such a non-adiabatic realization of 
quantum annealing. 

We next point out that a class of models of spin glasses has a special symmetry known as the 
gauge symmetry. We define the gauge transformation as the following local unitary operator for 
Pauli spin matrices, and the change of the sign of the interaction y,6j: 



Ga^G \G-X{G,, Gi-[ (| = -1) 



where a = x,y,z and Tij = Tij^i^j, where ^ is a classical gauge variable taking only ±1. 
After the above gauge transformation, the time-dependent Hamiltonian ([T]) is invariant, but the 
distribution function ^ is modified into P{Tij) = exp(/3pJTy^j^j)/2cosh/3pJ. This property 
helps us to derive the following results. 

3. Analysis 

3.1. Jarzynski equality 

Let us introduce another piece of theoretical tool, the Jarzynski equality, which is a generalization 
of the second law of thermodynamics. We here use the quantum version of the Jarzynski equality 
for the ±J Ising model, not the original version for the classical case [11^ I12j . with a specific 
configuration {rjj } as follows [T71 [18] , 

(e-^^)o.T = ff^, (4) 



where Z^(t;{rjj}) is the partition function for the instantaneous Hamiltonian, and /3 is 
the inverse temperature. The performed work during a nonequihbrium process is given by 
the difference between the outputs of the measurements of the initial and final energies as 
W = Em{T) — En{0), where m and n denote the indices of the eigenstates as H{T)\m{T)) = 
EmiT)\m{T)) and H{0)\n{0)) = £^„(0)|n(0)). The left-hand side of equation (gD expresses the 
average of the exponentiated work over all the realizations of nonequihbrium processes starting 
from the equilibrium ensemble, and is evaluated as 



(e-^^)o^r = 5:e-^(^'"(^)-^"(0))p(m|n) 

m,n 



m\n) = — - — 



where we use the path probability for a nonequihbrium process as, 
j(m(T)|[/o^T|f^(0))p. Here Uq^t is the time evolution operator. 



(5) 



P{m\n) 



3.2. Jarzynski equality for quantum annealing : Symmetric distribution 

The partition function for the initial Hamiltonian H{{)) does not depend on the configuration 
{rjj} since it is given only by the transverse field, Z^(0;{rjj}) = (2cosh/?r) . Here is the 
number of spins. Let us now consider the configurational average of the Jarzynski equality for 
quantum annealing by summation over all the possible configurations {Tjj} as 



o-s>r 



_ [Zp{T-{n,})]^^ 

(2cosh/3r)^ 

The explicit expression of the quantity on the right-hand side is 
[Z0{T;{Tij})]^^ I 



(2 cosh pvf (2 cosh /3r)^ (2 cosh l^pjf^ 



nj {ij) 



(6) 



(7) 



where Nb is the number of interactions. If we set /3p = 0, which corresponds to the case with 
the symmetric distribution of the interaction {p = 1/2), the above quantity can be evaluated as 



-I3W 



_ 2^ (2 cosh /3 J) 
/3p=o ~ (2cosh/3r)^' 



(8) 



3.3. Inverse statistics on special suhspace 

Let us show that the gauge transformation reveals a non-trivial property of the nonequihbrium 
average of the exponentiated work during quantum annealing. We take the configurational 
average of the inverse of the nonequilibrium-averaged exponentiated work as 



1 



(e-/3^)o^T 



_ (2 cosh I3T) 
/3p (2cosh/3pJ)^ 



N 



E 



rife-) 



Zp{T;{Tij})' 



(9) 



Here we apply the gauge transformation to the above equality. The quantity on the left-hand 
side does not change, since the Hamiltonian and the time-evolution operator are invariant. On 
the other hand, the quantity on the right-hand side becomes 



1 



O-i-T 



/3p 



(2 cosh I3T) 
(2 cosh /3p J) 



N 



Nb E 

Tij 



Zf)iT;{Tij}) 



(10) 



Setting 13 = (3p and taking the summation over all configurations of leads us to 



1 1 _ (2 cosh /3r) 



N 



(e-/5'^)o^TJ/3 2^ (2 cosh /3 J) 



Nb ' 



(11) 



which implies that 



-1 




1 



(12) 



This equality represents a non-trivial relationship in quantum annealing driven toward two 
different states in spin glasses /3p = and j3p = j3. It is important to remark that we do not 
use any approximations in the above calculation, and we deal with the nonequilibrium average 
during quantum annealing, which does not assume to be quasi-static, unlike quantum adiabatic 
computation. 

4. Summary 

We can obtain several further exact equalities for the repetition of non-adiabatic quantum 
annealing and establish non-trivial relations that hold between quantum annealing toward 
different spin glasses, as will be announced shortly in a separate paper. Exactness of our results 
would be useful as a benchmark to check experimental conditions of nonequilibrium process in 
quantum spin dynamics. Notice that our results are scalable, which means independence of the 
system size, since the Jarzynski equality holds for any size if we take all fluctuations into account 
in estimation of the expectation value on the left-hand side of equation ([!]). Finite errors on the 
ratio of the partition functions may be attributed to the rare event on nonequilibrium processes. 
The present study will help us to understand and manipulate quantum dynamics efficiently. 
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